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Abstract 

An upper limit on n — n oscillations can be obtained from the stability of matter. 
This relation has been worked out theoretically and together with data yields Tnn > 
(1 — 2) X 10*^ sec. A recent publication claims a different relation and finds a nuclear 
suppression of the n ^ n oscillations which is two orders of magnitude weaker than 
the previous evaluations. Using the same approach we find, nevertheless, that the 
earlier estimates are correct and conclude that future experiments with free neutrons 
from reactors are capable to put a stronger limit on the n — n oscillation time than 
experiments with large amount of neutrons bound in nuclei. 
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1. Introduction 

Although the proton decay expected in a simple version of grand unified theories has 
not been observed with a large lower limit on the lifetime, other variants of the theory 
predict reactions with nonconservation of baryon minus lepton number by two units, like 
neutron-antineutron oscillations or neutrinoless double-/? decay (see e.g. 0). The former 
phenomenon has been searched for with free neutrons from reactors in the recent experiment 
ILL-Grenoble and the lower limit for the oscillation time 

Tnn > 0.86 X 10^ sec (1) 

is established. A new experiment planned at ORNL is intended to increase this limit by 
two orders of magnitude. 

Experiments which search for the proton decay are sensitive ton ^ n transitions as well, 
since an antineutron produced inside a nucleus annihilates and blows the nucleus up. Such 
events are not observed above the background in the experiments ^. The experiments 
1^, P established quite high lower limits on the lifetime of the nuclei, respectively, 

T(0^^) > 4.3 X 10^^ (2) 

T(Fe^^) > 6.5 X lO^^yr (3) 

These data can be used to deduce a limit on Tnn provided the nuclear effects are known. 
Two groups and have calculated the nuclear effects and based on these data have 
provided close estimates of 

> (1-2) X lO^sec . (4) 

Recently, these calculations have been challenged and via the nonrelativistic diagram 
technique a value of Tnn is deduced from Eqs (0)-(||), which is about two orders of magnitude 
higher than Eq. (1). If the latter calculation were true, this lower limit would be comparable 
to what is planned for future experiments with free neutrons [Q. 

In this note we discuss this controversy and conclude that the theory which leads to 
@) is correct. We suggest an intuitive physical picture of a possible n ^ n transition in a 
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medium and point out why the approximation used in |^ is not vahd. In our calculations 
we emphasize that at low energy the annihilation radius is quite large which may lead to 
substantial corrections to (|). 



2. Interference effects in n — n oscillations 

We assume that the physical neutron (antineutron), i.e. the vacuum eigenstate, contains 
an admixture of a n (n) component. Suppose that at time t = we have a pure \n) state 
(neutrons from a reactor, see below), which is a mixture of the physical \n)pfiys and \n)phys 
having different masses. As a result of a phase shift between these states which increases 
with time a \n) component appears with a probability proportional to the square of time 
(compare with the well known effect of — oscillations). 



(provided that t <^ 1/e)- Here e = l/r„jj is related to the off diagonal matrix element of the 
vacuum Hamiltonian, and Tnn is the n ^ n oscillation time. 

The quadratic time dependence is a direct manifestation of quantum coherence, a linear 
dependence would be the result in classical physics. This fact is extremely important for an 
understanding of the suppression mechanism n n transitions in nuclei. Let us consider 
a set up where a neutron propagates a distance L to a counter which can measure the 
antineutron component. In another set up we have such counters each separated by a 
distance L/N . Each counter filters away the antineutron component. In the first case we 
have a probability ^ [tLjvf, in the second set up Wl^^ ^ N {t L/NvY ^ Wl/N, to 
detect the n.Here v is the velocity of the neutron. The second case is closely related to the 
absorption in a nucleus, where each bound nucleon serves as a detector for the n component. 

2.1 The mean lifetime of an antineutron in a medium 

The mean free time of propagation of the antineutron through nuclear matter (without 




(5) 
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annihilation) is tfree = ^/{cannVp), where aann is the n — N annihilation cross section and 
p is the density of nucleons. The smallness of tfree is one important source of suppression 
of n — s> n transitions in a medium. A neutron which propagates through a nuclear medium 
produces an antineutron which annihilates with a rate 

dWfi 2 2 /^N 

= (TannVp e t j^^^ , (6) 

In the limit of small velocity one can replace (Tann by the imaginary part of the n — N 
scattering length 0, 

47r 

Vdann = Im UnN (7) 

m 

The annihilation rate turns out to be inversely proportional the scattering length 

dt AirplmanN 

This result is quite different from that in where the annihilation rate is found to be 
proportional to Ima^^Tv Below we show that this might be true only in the limit of very 
small density p which is not the case for a nuclear medium. 



2.2 The coherence and formation times 

There is another source of the n n suppression. In the presence of an external field 
which acts differently on n and n, the transition n — >■ fi is forbidden by energy conservation. 
It may happen only virtually, with a lifetime of the n fluctuation in the neutron, 

tc = ^ , (9) 
AE ^ ^ 

where AE is the energy splitting between n and n. In order to be detected the n fluctuation 
has to interact (annihilate) during its lifetime. 

Note that at the far periphery of a nucleus the neutron is quasi-free, i.e. AE —>■ 0, and 
one may conclude from @ and (§) that the n fluctuation lifetime has no restriction. This is 
not true because the neutron cannot stay at the periphery longer than l/ur, where uj is the 
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oscillator frequency (the shell splitting). Thus, one must replace effectively AE AE + uj. 
This is equivalent to inclusion of the formation time effects 0]. 

We conclude that the effect of coherence and formation time becomes the main source 
of n — i> n suppression in the case of low density, 

dWfi _ 47r p Im anN „^ 
dt ~ m {AE + Loy ^ ^ 

The (unrealistic) limit of low density is recovered by calculations in [^]. There the energy 
splitting between n and n in nuclear medium is taken to be AE = 8 MeV, the mean binding 
energy of a bound nucleon. It is assumed that the antineutron appears momentarily and is 
unable to feel the mean field of the whole nucleus, therefore, it has no binding. This cannot 
be true because the n fluctuation lifetime (^) is very long in this case, = 25 fm. 

3. The evolution equation for a bound neutron 

To incorporate all the above effects one should consider the evolution of a n — n system 
in a medium. This is a typical two-channel problem (see e.g. in [0, 0]). 

tj^Mt,r)) = H{f)mt,r)) , (11) 

where the wave function |$(t,r)) contains the two components n and n and depends on 
time t and the space f. The evolution operator H{r) has a form. 



H(r) 



(12) 



y e Enir}-^p{r)lmanN J 
Here p(r) is the nuclear density modified due to annihilation radius Rann which may be 
large. 



The modified density p(r) spreads to larger distances than p(r). 
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The solution of equation ( pT]) for the weight of the antineutron component at long times 

(t > tctfree) reads, 



^) [A£;(r) + u;]2 + ((27r/m)p(r)Ima^Ar)2 ^^^^ 
We corrected here for the oscillator frequency considered in the previous section. 

One should average expression ([l^) over different neutrons which have different binding 
energies dependent on a shell [0. However, this difference is much smaller that the binding 
energy of an antineutron used in |]^, |^ which was assumed to be independent of the shell. 
Therefore, one can safely use the mean value of AE in (H) and then sum over the neutrons. 



The annihilation rate of antineutrons produced in a nucleus is given by integral of ( |I^ ) 
over the nuclear volume weighted by neutron density and annihilation cross section, 

(15) 



Tr AirlmanNJ p(f) 1 + K^{f 
Here A and N are the atomic and neutron numbers of the nucleus, respectively, and 

K(r-, = , ,^^<,^|" (16) 

(27r/m)p(r) |ImafiAr| 

The assumption made in ||^ corresponds to K{r) ^ 1. This could be not a very rough 
approximation, but AE = 8 MeV used in is an order of magnitude smaller than in |]^, |] 
and here (see discussion at the end of previous section). This is the main reason why Ta is 
overestimated by two orders of magnitude in 0. 

4. Numerical evaluation 

Following to [0] we assume that the n binding energy depends on r in the same way 
as p(r). Therefore, ([T6|) is approximately a constant K{r) ^ K{0) at r < R^, where 
Ra ~ tqA^/^ is the nuclear radius. The ratio p{r)/p{r) in (p!5|) is unity at r < Ra, but 
then steeply falls down. Using these observation we can perform a very simple (but rough) 
evaluation of (|T3p, 
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A spectacular observation which follows from this result is A-independence of the nu- 
clear annihilation width. Of course, better calculations (see below) may lead to a weak 
A-dependence. The claim in [0 that annihilation rate gets main contribution from the neu- 
trons on the surface might be misleading. In fact their result is also nearly A-independent, 
what supports our observation that all the nuclear volume contributes (see also in |^). 



The pp scattering length was determined in [O using the LEAR data on pp atoms, 
llmoppl f» 0.7— 1.2 fm. Following P we use the same value for nn and np. With this value 
we have, 

2n 

— p(0) llmanNl ^ 30 - 50 MeV , (18) 
m 

where we use the central density p(0) = 0.16 fm~^. 

Expression (|T^) corresponds to the imaginary part of the antineutron nuclear potential 
in 1^. The real part of the antineutron potential in the center of the nucleus can be only 
guessed since the data on antineutron atoms are sensitive only to the surface of the nucleus. 
This is the reason of diversity of antineutron potentials used in [0, which we consider 
as very unreliable. An advantage of and our approach is that we use experimental 
information on nn scattering length, rather than guess a value of imaginary part of the 
nuclear potential. For a simple evaluation we choose the real part to be of the same order 
as the imaginary one, i.e. that K{0) ~ 1 (A£^(0) ~ 40 MeV). This leads to the estimate 
for the nuclear factor Tr in (|T^), Tr ^ 0.23 /m~^, which is quite close to the results of 
exact integration in (|T3]): Tr = 0.20 /m^^ for ^^O and Tr = 0.18 /m~^ for ^^Fe. For these 



calculations we used the oscillator parameter lu from [0 and AO MeV for ([T8|). Following 
0, 1^ we assumed that AE{r) oc p(r). We use the realistic Woods-Saxon parameterization 
for the nuclear density. Dependence of the result on the annihilation radius turns out to be 
quite weak. We use Rann = 1 f^- 

We take the upper bound for AE{0) to be 200 MeV which exceeds most of model values 



used in 0, §]. Corresponding values of Tr are 1.2 fm~^ for ^^O and 1.3 fm~^ for ^^Fe. 

This estimation for nuclear effects combined with the experimental limits for the lifetime 
of the nuclei leads to the following restrictions for the oscillation time Tnn = 



These results are very close to the low bound (||) found in experiment with free neutrons 
and are compatible with (Q) - as well. 

The uncertainty of the lower limit (^) is still conventional since it depends on the 
guessed real part of the antineutron potential inside the nucleus. 

Concluding, we do not see how the lower limit for Tnn can be pushed up significantly 
by using the stability of matter and one may have to turn again to experiments with free 
neutrons with long fiight path and small velocity. 
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